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ABSTRACT 

We present a new mass estimate for th e Hercules dwarf s pheroidal galaxy (dSph), based on 
the revised velocity dispersion obtained by lAden et al 



The removal of a significant foreground contamination using newly acquired Stromgren pho- 
tometry has resulted in a reduced velocity dispersion. Using this new velocity dispersion of 
3.72 ± 0.91 kms , we find a mass of M 300 = l-Vtl.l x 10 6 M Q within the central 300 pc, which 
is also the half-light radius, and a mass of M433 = 3.7^'g x 1O 6 M0 within the reach of our data 
to 433 pc, significantly lower than previous estimates. We derive an overall mass-to-light ratio of 
M^/L — 103^4g [Mq/Lq]. Our mass estimate calls into question recent claims of a common 
mass scale for dSph galaxies. 

Additionally, we find tentative evidence for a velocity gradient in our kinematic data of 16 ± 
Skms^kpc^ 1 , and evidence of an asymmetric extension in the light distribution at ~ 0.5 kpc. 
We explore the possibility that these features are due to tidal interactions with the Milky Way. 
We show that there is a self-consistent model in which Hercules has an assumed tidal radius 
of r t = 485 pc, an orbital pericentre of r p = 18.5 ± 5 kpc, and a mass within r t of Aftid,r t = 
5.2^2 7 x 10 6 M Q . Proper motions are required to test this model. Although we cannot exclude 



1. Introduction 

Dwarf spheroidal (dSph) galaxies are believed 
to play an important role in the formation and 
evolution of much more luminous galaxies (e.g. 
Gallagher &: Wvsd I1994 ). dSphs are character- 
ized by their low surface brightness, low total lu- 
minosity and spheroidal shapes which are consis- 
tent wit h their pressure-sup ported stellar kinemat- 
ics (e.g. 



Grebel et all 12003) 



Knowledge of dSph masses is essential for com- 
parison with cosmological simulations of galaxy 
formation. Good mass estimates help us to es- 
tablish whether the paucity in the number of 
observed systems (a few tens) to the number of 
predicted satellite haloes (several thousands) rep- 
resents a fun damental fa i lure o f our cosmological 
model (e.g. Moore et a l. 1999), or whether it is 
simply telling us that galax y formation is inef - 
ficient on small scales (e.g. iRead et aL 2006a). 
Recent studies have suggested that the dSph 
galaxies share a common mass within a certain 
radius (IWalker et all l2007t IStrigari et all 120081 : 
Walke r et al. I l2009h . If confirmed, this would be 
an important clue to the processes which regulate 
the formation of the lowest luminosity galaxies. 

All mass estimates implicitly assume that con- 
tamination of the kinematic sample by foreground 
stars or unbound tidal tails is negligible, and 
that the system is in (or close to) virial equilib- 
rium. Since the mass of an equilibrium stellar 
system is proportional to its velocity dispersion 
squared, an over-estimate of the velocity disper- 
sion will result in an inflated mass estimate. The 
assumption of virial equilibrium may be called 
i nto question if the system i s tidally disrupting 
(lOh et al.lll995HKroupa|l997tlKTimentowski et al 
20071: iMunoz et al.ll200; 



foreground dwarf stars. By weeding out the fore- 
ground contaminants we found that the dispersion 
for Hercules is reduced from 7.33 ± 1.08 kins -1 to 
3.72 ± 0.91 kms -1 . In this letter we explore the 
consequences of this finding. 

In Section 2, we derive a new mass for Hercules 
and show that it is not consistent with a common 
mass scale for the dSphs. In Section 3, we in- 
vestigate the possibility of a velocity gradient in 
the kinematic data of Hercules. In Section 4, we 
discuss the relative importance of tides and dark 
matter in Hercules. Section 5 summarises our con- 
clusions. 

2. Mass estimate from the spherical Jeans 
equations 

We use t he velocity disper sion, a v = 3.72 ± 



The recently discovered Hercules dSph dBelokurov et al 



2007), with an cllipticity of e ~ 0.5 (| Coleman et al 



2007[), is an example of a galaxy in which as- 
sumptions o f equ ilibrium may be incorrect. In 
Aden et al. ( 20091 ) we showed that the mean ve- 
locity of the Hercules dSph is embedded in the 
foreground dwarf star velocity distribution. In 
order to weed out the foreground dwarf stars 
we used the Stromgren c\ index. This index is 
able to clearly identify red-giant branch (RGB) 
stars redder than the horizontal branch, enabling 
a separation of RGB stars in the dSph galaxy and 



0.91kms _1 (jAden et a l. 2009), to constrain the 



mass of the Hercules dSph. We assume that the 
system is in dynamical equilibrium, is spherically 
symmetric, has an isotropic velocity distribution, 
and a flat velocity dispersion profileo With these 
assumptions, the Jeans equ ation for the mass 
distri bution (Eq. (4.215) in Binnev fc Tremaind 



20081 ) becomes: 



dv{r) _ v(r)GM{r) 
dr r 2 



(1) 



where r is the three-dimensional distance from the 
centre of the galaxy, v{r) is the de-projected stel- 
lar density profile and M(r) is the e nclosed mass. 
We us e the exponential profile from iMartin et aT 
(2008) to describe the stellar density profile 



The de-projected exponen tial profile is given by 



Klimentowski et alj (|2007l ) as (setting to = 1 in 



their Eq. (4)): 



v(r) = v a 



-0.445 



-r I a. 



(2) 



where a is the exponential scale radius, related to 
the half-light radius as r/j = 1.68 a, and vq is the 
central stellar density. Solving Eq. |T]) for M(r) 
using Eq. ([2]) yields 



M(r) 



r(r + 0.445 a)a% 



(3) 



!Note that IWolf et al.l (I2009T ) and IWalker et alj (I2009D find 
that their mass estimates at the half light radius are insensi- 
tive to a wide range of mass models and velocity anisotropy 
parameterisations, so our results should not be sensitive to 
these assumptions. 
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Using a half-light radius of 330I52 P c <|Martin et al 
2008T ) we solve Eq. ((3|) for r = 433 pc, which corre- 



sponds to the outermost member in our kinematic 
sample (Tabled]). 

We estimate the error in the mass using 10 5 
Monte Carlo realisations of the half-light radius, 
distance and velocity dispersion of Hercules drawn 
from within the individual error bars on each pa- 
rameter. In this way, we obtain M433 = 3.7^'g x 
10 6 Mq (Fig. [T|). The quoted errors are 1-a lim- 
its from the Monte Carlo sampling. Assuming 
a total luminosity of L = (3.6 ± 1.1) x 10 4 L r . 



( Martin et al. 20081 ). we find a median mass-to- 



light ratio of M 433 /L = 103l|| [Mq/ L 



Aden et al. ( 20091 ) emphasise the importance 
of weeding out foreground, contaminating dwarf 
stars. This is particularly vital for the Hercules 
dSph as its systemic velocity coincides with the 
bulk motion of dwarf stars in the Milky Way disk 
in the direction of Hercules. Using the contami- 
nated^! velocity dispersion (Table []} we would ob- 



+0.5 
-0.4 



IO'Mq, almost a factor 



tain a mass of 1.4 
of four larger than the uncontaminated estimate. 

2.1. Hercules and the "Common Mass 
Scale" 



Strigari et al.1 (|2008f ) speculate that the mass 
within a fixed radius, 300 pc, is approximately the 
same (~ 10 7 Mq) in all the observed dSphs. Our 
revised velocity dispersion implies a mass within 
the inner 300 pc (M 300 ) of only 1.9±oi x 10 6 M Q . 
This indicates that Hercules falls considerably be- 
low this "common mass scale" for dSphs. Since 
300 pc is also approximately the half-light radius, 
this implies that Hercule s also lies significantly off 
the IWalker et all (|2009l ) enclosed half-light mass 
scaling relation. 

To confirm that our result is not sensitive to 
our choice of surface brightness profile, we have 
repeated our calculation using a Plummer pro- 
file with scale rad ius r p i — 321pc (as used by 
Strigari et al. 20081 ). In this case we find M 30 o = 
2.3 x 10 6 Mq in agreement with the mass cal- 
culated above. If instead we use our contami- 



2 Note that this "contaminated" dispersion has had its 3er 
velocity outliers removed. Also, the initial candidate selec- 
tion for the velocities were chosen using colour-magnitude 
diagram cuts. Even so, a significant fraction of foreground 
stars remain in this sample and are detected only through 
the use of the Stromgren photometry. 



nated velocity dispersion and the exponential sur- 
face density profile, we obtain M 30 o = 7.4+ 2 'J x 
10 6 M Pl which agree s with the value presented in 
Strigari et all (|2008h . 

However, we note that IStrigari et al. (2008) 
used a dispersion of 5.1 ± 0.9 km s -1 (taken from 
Simon fc Gehal 120071 ) to obtain their mass esti- 
mate. This is smaller than our contaminated value 
of 7.33 ± 1.08kms _1 , yet they obtain a similar 
median mass to us. If we use a dispersion of 
5.1 ± 0.9kms _1 as they do, we obtain a mass 
M 300 = 3.6ti 2 x 10 6 M o . This is just consistent 
with their determination within our mutual error 
bars. We note that the median M300 will depend 
on the details of the mass modelling procedure and 
so can be expected to differ between their study 
and ours (Matt Walker priv. comm.). 

Additionally, we have repeated our calculation 
of the mass within 300 pc using a more recent es- 
timate of t he half-light radius, — 230 pc, by 
Sand et al. (2009). This half-light radius is smaller 
and gives a mass that is ~ 1.5 times more massive 
than the mass calculat e d usin g the half-light ra- 
dius from [Martin et alj (|2008l) . 



3. A velocity gradient in Hercules 

The presence of a velocity gradient in a dSph 
could either be indicative of an intrinsic rotation, 
or a sign of tidal interaction with the Milky Way. 
In this section, we test for possible velocity gradi- 
ents in Hercules. 

Assuming that the rotation around the semi- 
minor axis (in the ellipse that desribes the ori- 
entation of the dSph) is more likely then around 
the semi-major axis we can define the "semi-minor 
axis distance" d m i(9) as the distance, perpendicu- 
lar to the semi-minor axis, between the axis and 
the star. For each of the 18 RGB stars with radial 
velocities we calculate this distance for different 
position angles 6. We then derive the gradient , 
fc rot , and zeropoint, m„, by solving the following 
equation with a least-squares fit for each position 
angle, 

Vrad = fc rot X d mi (6>) + m v (4) 

where V, a d is the radial velocity measurement for 
each stars. The least-squares fit to this function 
yields a x 2 value for each (Fig. [2k). 

The distribution enclosed by Xm ; - + 1 corre- 



sponds to 1 a for a normal distribution ( Press et al 
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1992). We use this to obtain the error in the po- 
sition angle which minimizes x 2 ■ 



We find a position angle of — 35I 2 3 degrees with 
a velocity gradient of 16 ± 3kms _1 kpc _1 , and a 
zeropoint of 45.11 ± 0.38 kins" 1 . We obtain a re- 
duced x 2 of 3.89 for our 18 stars with 2 degrees of 
freedom. 

( 2008h . 



Following IWalker et al 



we 



estimate 

the significance of the velocity gradient using 10 5 
Monte Carlo realisations. In each realisation, we 
sample the velocity and spatial distributions in- 
dependently. Thus we scramble the correlation 
between velocity and position, while maintaining 
the original velocity distribution and spatial posi- 
tions. This is valid as long as the phase space dis- 
tribution function of the stars is separable (which 
implies that the velocity dispersion is independent 
of radius). We determine the significance of the 
velocity gradient by computing the fraction of re- 
alisations that fail to produce a \ 2 as l° w as the 
one calculated from the real data. We find a sig- 
nificance of the velocity gradient of 78 per cent 
(1.23c-). 

4. Discussion 

4.1. Galactic tides in Hercules? 

In the previous section we found tentative evi- 
dence for a velocity gradient in the Hercules dSph. 
The peak-to-pcak difference of ~ 10 kins" 1 within 
a radius of less than 1 kpc could be interpret ed as 
the effect o f Galactic tides ( Read et al ]|2006bh . In- 
terestingly, iMartin I (|2009h recently estimated the 
orbit of Hercules based on the assumption that its 
elongation is tidally induced, and predicted a ve- 
locity gradient of at least 7kms" 1 kpc~ 1 which is 
consist ent with our obser ved gradient. Addition- 
ally, in lAden et all (|2009h we found that the spa- 
tial distribution of the Hercules stars is asymmet- 
ric at ~0.5 kpc. There are three significant outlier 
stars to the South, but no corresponding members 
at this distance in either the Northern or Western 
fields (Fig. [3]). These three stars are unambigu- 
ously identified as RGB stars from Stromgren pho- 
tometry. We now consider the possibility that the 
velocity gradient and the positional outliers are 
evidence that Hercules is being tidally disrupted, 
and use this information to obtain a second mass 
estimate for the system. 

The tidal radius of a dSph depends on the 



potential of the host galaxy, the potential of 
the dSph, the orbit of the dSph within the host 
galax y and the orbit of the stars within the dSph 



(e.g. iRead et al.l l2006dl . We parameterise the 



Milky Way pot e ntial using the default model in 
Johnston et alj (|2005i ). analysed in the Galactic 
plane; and the Hercules potential using a gener- 
alised Hernquist profile ( Hernquistlll99C 



7—4 



(5) 



where M is the total mass, r s the scalelength, and 
7 the central logarithmic cusp slope. We consider 
the ranges 0.3 < r s < 3 kpc and < 7 < 1. Our 
results are not sensitive to these choices. 

The orbit of the Hercules dSph is unknown 
but it currently lies at a distance of 132 kpc from 
the galactic centre, and has a heliocentric veloc- 
ity of 45.2 ± 1.09 km s" 1 , which implies a galac- 
tocentric r adial velocity of 145 ± 1.09 kms" 1 (see 
Eq. (5) in lCourteau fc van den Berghlll999l ). Us- 
ing the above potential model for the Milky Way 
and setting the tangential velocity component for 
Hercules to zero, this gives us a minimum apoc- 
entre for Hercules of r a , m i n = 188.5 kpc. Thus, 
we consider apocentre and pericentre ranges of 
188.5 < r a < 600 kpc and 5 < r p < 132 kpc, re- 
spectively. 

IRead et"ahl ( 2006d ) find that photometric fea- 
tures are typically seen beyond the retrograde 
tidal radius. To proceed, we make the assump- 
tion that the outliers in Fig. [3] indicate the lo- 
cation of the retrograde tidal radius of Hercules, 
i.e. r t ~ 485 pc. If we assume further that the 
tidal radius of Hercules is set at th e pericentre 
of its o rbit, we can solve Eq. (7) of IRead et al. 
(|2006dl to calculate the mass, M t id, of Hercules as 



a function of its orbital pericentre r p , for our as- 
sumed ranges of r s , 7 and r a . We calculate both its 
mass within r t (M t id,r t ) and its mass within 433 pc 
(M t id. 433) which can then be compared with our 
mass estimate based on the spherical Jeans equa- 
tion (Msj,433). The results of this calculation are 
given in Fig. 0J The horizontal light grey band 
marks Msj,433, the dark grey band marks Mtid, 433 
and the medium grey band marks M t id,r t ■ The 
vertical solid and dashed lines mark the pericen- 
tres at which M t id,433 = Msj,433. The width of 
the tidal mass bands is due to the unknowns: r s , 
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7, and r a . However, to lowest order rt depends 
only on the mean density enclosed within it, and 
so these bands are narrow. For this reason, we ob- 
tain an estimate of both the orbital pericentre of 
Hercules (r p = 18.5 ± 5kpc) and its mass within 



the tidal radius (Aftid,r t = 5.2^_ 2 7 



h2 - 7 x 1O 6 M ). The 
primary source of error on both of these quanti- 
ties is our assumed tidal radius r t . Empirically we 
derive scalings of: 



32 



n 



1 kpc 



kpc 



(6) 



with (i — 0.76 over the range r t — [0.4, 2] kpc, and: 



tid,r t 



= 24 x 10 b 



(7) 



1 kpc 

with a = 2.1 over the same range. 

If we assume rt = 485 pc, we obtain a similar es- 
ti mate of the peri centric distance as that obtained 
bv lMartinl (|2009l ). 



4.2. Is Hercules a dark matter free sys- 
tem? 

We have shown that our Hercules data are con- 
sistent with the presence of dark matter. We 
now consider whether models without dark matter 
could also reproduce the data. 

The most extreme scenario is that Hercules is 
disintegrating and its velocity dispersion arises 
solely due to the motion of its unbound mem- 
ber stars. If so, it will rapidly become unobserv- 
able (i.e. reach a lower surface brightness than 
the detection limit of the SDSS). Hercules has 
a surface brightness of 27.2 ± 0.6 mag arcsec -2 



( Martin et al. 20081) . If the unbound stars are 



moving away from Hercules at a velocity equal 
to the velocity dispersion (3.72 ± 0.91 kms -1 ) it 
would require only ~ 200 x 10 6 years for the Her- 
cules dSph to fall below the detection l i mit o f 

2008). 



Koposov et al 



SDSS (~ 30 mag arcsec" 2 
Given this short timescale, it is very unlikely that 
we would observe Hercules at this phase of its evo- 
lution. 

Fellhauer et al. (2007) simulated the disruption 
of the UMa II dSph. Using a model that does 
not distinguish between luminous and dark mat- 
ter they simulate the surface brightness, velocity 
dispersion and the mean radial velocity for UMa 
II after 9, 10 and HGyr (their Fig. 9). 



In the absence of dark matter, the equilibrium 
velocity dispersion (i.e. a 2 — G M/r) for the Her- 
cules dSph would be ~ lkms -1 , assuming a stel- 
lar mass of 5 x 10 4 M Q and a radiu s equal to the 



obser ved half-light radius (300 pc; iMartin et al 
l2008h . This is much lower than the measured ve- 
locity dispersion of 3.72 ± 0.91 kms -1 . We con- 
clude that the Hercules dSph can only be a dark 
matter free system if its velocity dispersion has 
been inflated significantly - i.e. if it is in the ad- 
vanced stages of tidal disruption. 



However, in the simulations of Fell hauer et al 



(|2007l ) the evolutionary phase in which a tidal rem- 
nant exhibits both an inflated velocity dispersion 
and a velocity gradient, and remains centrally con- 
centrated is very short. Within ~ 1 Gyr the sys- 
tem goes from bound, with almost no signs of tidal 
disturbance, to complete disruption. 



(We note that iKroupal (|1997I ) performed sim- 
ulations of dSph galaxies without dark matter, 
and found that it is possible to obtain long-lived 
remnants whose properties are remarkably similar 
to those of Hercules in terms of velocity disper- 
sion and velocity gradient (see his Fig. 9) . How- 
ever, the remnants simulated by IKroupal (|l997l) 
are significantly more luminous th an the Hercules 
dSph. Based on the simulations of lFellhauer et all 
(l2007h . it seems unlikely that lower luminosity, 
purely stellar remnants with the correct proper- 
ties could survive for a significant time.) 

It seems difficult to understand Hercules with- 
out dark matter. However, it could be dark mat- 
ter dominated and experiencing significant tidal 
disturbance. In this case, our simple Jeans analy- 
sis may have over-estimated its mass and Hercules 
might lie even further away from a possible com- 
mon mass scale for the dSph galaxies. 

5. Conclusions 



We have calculated the mass of the Hercules 
dSph using the n ew velocity dispers ion for the sys- 
tem obtained by I Aden etail (|2009h . We find that 
the mass within the volume enclosed by our ob- 
served stars is 3.7^ 2 'g x 10 6 M , leading to a mass- 
to-light ratio of 103^ [Mq/Lq]. Interestingly, 
the mass within 300 pc is significan tly lower than 
the "c ommon mass scale" found by Strigari et alJ 
([20081 ). suggesting that Hercules does not share 
the halo properties seen in other dSphs. 
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We found tentative evidence for a velocity gra- 
dient of 16 ± 3kms _1 kpc" 1 , and evidence of an 
asymmetric extension in the light distribution at 
~ 0.5 kpc. We explored the hypothesis that these 
features are due to tidal interactions with the 
Milky Way. Assuming a tidal radius of 485 pc, we 
show that a self-consistent model requires Hercules 
to be on an orbit with pericentre r p — 18.5 ±5 kpc, 
and with a mass within r t of A/ t id,r t = 5. 2^2 ^ x 
10 6 M o . 
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ing discussions. D.A thanks Lennart Lindegren at 
Lund Observatory for help with statistics. S.F. is a 
Royal Swedish Academy of Sciences Research Fel- 
low supported by a grant from the Knut and Alice 
Wallenberg Foundation. M.I.W is supported by a 
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Fig. 1. — Distribution of Monte Carlo mass esti- 
mates for Hercules. Solid line indicates the median 
from the Monte Carlo sampling. Dashed lines indi- 
cate the la limits from the Monte Carlo sampling, 
(a) Mass within the radius defined by the outer- 
most RGB star in the kinematic sample, r — 433 
pc. (b) Mass within the inner 300 pc. 
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Fig. 2. — (a) x 2 vs - the position angle, (b) V ra d 
vs. semi-minor axis distance. The solid line in- 
dicates the least-square fit. • indicates members 



based on Stromgren photometry (see lAden et al 
2009h . 



Fig. 3. — Spatial distribution on the sky of the 
stars identified as Hercules members from both 
spectroscopy and photometry, o indicate mem- 
bers based on Stromgren photometry. • indi- 
cate members based on Stromgren photometry for 
which there are radial velocity measurements. The 
dashed line indicates the semi-minor axis for the 
position angle of our detected rotation axis. The 
solid ellipse represents the core radius and the dot- 
ted ellipse the in ner border of the field region of the 



Hercules dSph ([Coleman et alJ 120071 ). Solid lines 



outline the footprint of the WFC used to obtain 
the Stromgren photometry. 
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Fig. 4. — Tidally determined mass for Hercules 
as a function of orbital pericentre r p . We as- 
sume a retrograde tidal radius of rt = 1 kpc for 
this calculation. The horizontal grey band marks 
our mass estimate from the spherical Jeans equa- 
tion Mgj,433 (See Sect. 2). The dark grey band 
marks our tidally determined mass within 433 pc 
Aftid.433 for Hercules. The medium grey band 
shows Mtid,r tJ our tidally determined mass within 
r t . The vertical solid and dashed lines mark the 
pericentres at which M t id,433 = Afgj,433 (see Eq. 
[5] for a formula giving the scaling of M ti( j n with 
rt)- 
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Table 1 

Summary of determination of systemic velocities, velocity dispersions, masses and 

metallicities for the hercules dsph. 





Number of stars 


Vsys 

[kms^ 1 ] 


a 

[kms^ 1 ] 


Mpurthest 

[M@] 


Star 


M300PC 

[M S ] 


<[Fe/H]> * 
dcx 


Aden et al. (2009) 
This study 
This study 


28 [RGB stars, Ci sol.] 

32 [Only V rad sol.] 

18 [RGB stars with Vrad] 


40.87 ± 1.42 
45.20 ± 1.09 


7.33 ± 1.08 
3.72 ± 0.91 


1.4+3:5 x 
3.7tU x 


10 7 
10° 


7.4t 2 2 l x 10 6 
l-9±o,8 x 10 6 


-2.35 ± 0.31 



a See lAden et all J2009D 



Note. — Columns: (1,2) number of stars in each study and a short description of how the stars were selected; (3,4) systemic 
velocities and velocity dispersions; (5) mass within the ra dius defined by the outermost RGB star in our study; (6) mass within the 
inner 300 pc; (7) mean metallicity [sec Aden et al. (2009) for a full discussion of metallicities] . 



